Robust synchronization of an arbitrary number of spin-torque driven vortex nanooscillators 



Sergey Erokhin and Dmitry BerkovQ 

Innovent Technology Development, Priissingstrafie 27B, D-07745, Jena, Germany 
(Dated: February 8, 2013) 

Non-linear magnetization dynamics in ferromagnetic nanoelements excited by the spin-polarized dc-current 
is one of the most intensively studied phenomena in solid state magnetism. Despite immense efforts, synchro- 
nization of oscillations induced in several such nanoelements (spin-torque driven nanooscillators, or STNO) still 
represents a major challenge both from the fundamental and technological point of view. Here we propose a 
system where synchronization of any number of STNOs in a wide current range can be easily achieved, and 
the inclusion of thermal fluctuations does not disturb it. Our STNOs are represented by vortices inside squared 
nanoelements, which precession is driven by current injected via point contacts. Using full-scale micromagnetic 
simulations we show that synchronization of these STNOs is extremely dynamically stable due to a very large 
coupling energy provided by the magnetodipolar and exchange interaction between nanoelements. Finally, us- 
ing as an example a system of 10 oscillators, we demonstrate that our concept allows robust synchronization of 
arbitrary number of STNOs, even when corresponding nanocontacts have a broad size distribution. 



Magnetization excitation and switching in thin ferromag- 
netic films and nanoelements induced by the injection of a 
spin-polarized dc-current is a very active research area of 
modern solid state physics due to the large importance of these 
phenomena for our fundamental understanding of electrical 
and magnetic properties of condensed matter [1-3]. Particu- 
larly, an enormous attention has been paid to excitations of a 
special type - steady state gyration of magnetic vortices gov- 
erned by the spin-polarization of injected current. This ob- 
ject is a highly non-trivial example of a nanosized oscillator 
with strongly nonlinear properties. For this reason the study 
of its dynamical behavior, which should lead to the decisive 
improvement of our general understanding of interaction ef- 
fects in such systems, is one of the most important tasks of 
the current research in non-linear magnetic phenomena. Espe- 
cially various synchronization regimes in ensembles of these 
oscillators represent a very challenging fundamental problem, 
which solution for various cases has been the subject of many 
recent publications Q4J-42D - 

This interest is strongly supported by many promising ap- 
plications of this novel effect fl 1 0L lllll. Among these appli- 
cations, one of the most important would be a nanosized mi- 
crowave signal generator which frequency could be tuned by 
changing either the external magnetic field or the dc-current 
strength. Such generator could be very useful in, e.g., detec- 
tors of microwave fields and telecommunication appliances. 

Despite large experimental and theoretical efforts aimed 
to create such generators, the power and linewidth of a sin- 
gle spin torque-driven nanooscillator (STNO) are still not 
competitive with corresponding non-magnetic state-of-the-art 
devices. Decisive improvement of STNO-based generators 
is expected from synchronization (phase-locking) of several 
STNOs, what would result both in a quadratic increase of the 
output power (~ N 2 for N synchronized oscillators) and sub- 
stantial narrowing of the generated linewidth. Unfortunately, 
although intensive theoretical research led a substantial deep- 
ening of our understanding of synchronization process of two 
and more STNOs iH-H LL2l [3, experimental achievements 
here remain rather modest: in pioneering reports Q14| |l5| two 



closely placed point contact STNOs could be synchronized, 
and afterwards only in one experimental paper [|16ll the syn- 
chronization of four such oscillators was demonstrated. 

In most theoretical and experimental studies synchroniza- 
tion of STNOs with the nearly collinear magnetization con- 
figuration under the contact is considered. Recently it has 
been shown that the steady- state gyration of a magnetic vortex 
11 71 Il8l1 is also a possible candidate for realizing a system of 



synchronized STNOs §00. 

In this study we present the system design which allows 
the robust synchronization of an arbitrary number of vortex- 
based STNOs in a wide range of current strengths, even when 
the diameters of corresponding point contacts are significantly 
different. The proposed system exhibits a very high coupling 
energy between the constituent STNOs. For this reason the 
synchronization is stable with respect to thermal fluctuations, 
as we have also verified by corresponding simulations. 

To provide a basis for the understanding of our many- 
vortices system, we start with the brief description of the mag- 
netization dynamics of a single vortex-based STNO. The vor- 
tex nanooscillator used in this study is based on the square- 
shaped nanoelement with lateral sizes 500 x 500 nm 2 and 
thickness h = 20 nm; this thickness corresponds to a typi- 
cal experimental system (see, e.g., lfl9ll ). The current with 
the out-of-plane spin polarization direction and polarization 
degree P = 0.3 is injected via the point contact with the di- 
ameter D = 100 nm (unless stated otherwise, see last sec- 
tion) in the center of the nanosquare. Magnetic parameters 
typical for Permalloy are chosen: magnetization Ms = 860 
G, exchange constant A = 1.0 x 10" 6 erg/cm, negligibly 
small crystallographic magnetic anisotropy and Gilbert damp- 
ing X - 0.01. To prevent the reversal of the core polarity for 
high current values (see below for details), the constant ex- 
ternal field //perp = 2 kOe was applied perpendicularly to the 
nanoelement plane. This sample and all other systems stud- 
ied below were discretized only in-plane using a 4 x 4 nm 2 
mesh. Micromagnetic determination of the equilibrium sys- 
tem state and dynamic simulations of the spin torque-induced 
magnetization dynamics were carried out with our MicroMa- 
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FIG. 1. (color online) Magnetization dynamics of a single STNO: 
(a) and (b) - snapshots of typical magnetization configurations in 
low- and high-current oscillation regimes (red solid circles define 
the nanocontact area, the green dotted line in (a) is the vortex tra- 
jectory); (c) spectral amplitude of magnetization oscillations under 
the contact. 



gus software package 112011 , which dynamical part solves the 
Landau-Lifshitz- Gilbert equation of motion for magnetic mo- 
ments including the additional term for the spin- torque. 

For all systems studied in this article below magnetization 
dynamics was simulated during 200 ns of the 'physical time', 
starting from the time moment when the current was switched 
on. In order to explore only the steady precession regime, all 
spectra presented below have been calculated from the data 
collected for 100 ns < t < 200 ns. 

The next important issue is how to initiate the vortex gyra- 
tion, because without the external in-plane field and due the 
presence of the circularly symmetric Oersted field, the vortex 
core position in the center of the nanocontact corresponds to 
the state with the minimal magnetic free energy. To induce the 
vortex oscillations, an in-plane field pulse with the amplitude 
200 Oe and duration 3 ns was applied. This pulse was strong 
enough to 'knock out' the vortex core from the nanocontact 
center to the region outside the contact. In this region the com- 
peting influences of the spin torque, Gilbert damping, Oersted 
field and the stray field from the borders of the nanosquare led 
to the core motion around the point contact. After the system 
'forgets' the initial in-plane field pulse, this rotation occurs by 
a nearly circular orbit. 

Simulation results for the single nanooscillator described 
above are shown in Fig. \Hc). Note, that on this and all other 
figures we display the spectral amplitude of oscillations of the 
in-plane (m x ) magnetization component under the point con- 
tacts), thus simulating the effect measured by GMR. Simu- 
lations performed in the current interval 0.5 < / < 10 mA 
reveal that there exist two dynamical regimes: 'simple' vor- 
tex rotation around the nanocontact (for / < 8 mA) and the 
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FIG. 2. (color online) Magnetization dynamics of two dipolar- 
coupled STNOs: (a) - snapshot of a typical magnetization configura- 
tion in a synchronized regime (note opposite oscillation phases under 
the two contacts); (b) time dependencies of the average magnetiza- 
tion under the 1st (m\(t), red line) and 2nd (m x 2 (i), blue line) contacts 
and the their sum (mf ot (t), green line); (c) - spectral amplitude of the 
sum of magnetizations under the contacts m* ot (t) = m x (t) + m x (t); 
(d) total oscillation power of the difference rn x m (t) = m x (t) - m x 2 {i), 
demonstrating the large power increase in the synchronized regime if 
the additional phase shift of n is introduced between the ac-currents 
generated by the two contacts. 



regime, where the rotating vortex generates vortex- antivortex 
(V-AV) pairs (/ > 8 mA). In the latter case, the V-AV pair 
annihilates very quickly again 12111 . leading to the emission 
of a spin- wave burst, which is the dominating mechanism of 
the energy dissipation in the high-current regime. For the spin 
torque-induced vortex dynamics, this process has been dis- 
cussed in detail for a somewhat different point-contact system 
in 11911 . so that we do not elaborate this issue further. We only 
note, that in our case there is no core polarity switching due 
to the out-of-plane external field, so that the remaining vortex 
simply increases its angular velocity again until the next V-AV 
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FIG. 3. (color online) Magnetization dynamics of the two STNOs 
connected by a 'bridge' (partial exchange coupling), shown in same 
way as in Fig. 0c). STNOs are synchronized almost in the whole 
current region, where their steady- state gyration is observed. 
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FIG. 4. (color online) Magnetization dynamics of the two STNOs 
with a full exchange coupling, presented in same way as in Fig. 0c). 
Vortex gyration is out-of-phase almost for current values, leading to 
in-phase magnetization oscillations under the contacts. 



pair is generated. We present a snap-shot of the m x -component 
for the first regime in Fig. Qla) (white corresponds to m x = +1 
and black - to m x - -1) and a snap-shot of the out-of-plane 
m-component immediately after the annihilation of the V-AV 
pair for the second regime - in Fig. 02b); in this last image, 
diverging spin waves generated by the V-AV annihilation can 
be clearly seen. 

Now we can proceed with the analysis of synchronization 
in a systems of two nanooscillators. For the synchronization 
of STNOs, the first (low-current) regime is obviously more 
suitable. We point out, that in this regime for our system 
the vortex rotation frequency depends on the current strength 
very weakly already for small currents and is nearly constant 
for currents above « 4 mA (see Fig. \Hc)). This behavior 
is in a strong contrast with the results reported previously for 
the in-plane magnetized nanocontacts 1141 ll5|] and also differs 
considerably from the /(^-dependence found in InlEal]. De- 
tailed discussion of possible reasons for this difference will be 
presented elsewhere. Here we would like only to emphasize 
that such a weak dependence /(/) is a good prerequisite for 
the synchronization of STNOs in a broad current range. 

In this study, we investigate three kinds of coupled STNO 
systems, where the oscillators are coupled in different ways: 
(i) dipolar coupling only, i.e., nanosquares (placed in one 
plane) are fully separated by a narrow gap having the width 
d = 48 nm, see Fig. Eta); (ii) dipolar and partial exchange 



coupling, i.e., nanosquares are partially connected by the 
bridge with the same thickness h = 20 nm and the width 
/ = 200 nm, made of the same magnetic material (Fig. E^a)); 
and (iii) dipolar and full exchange coupling, realized as a sys- 
tem of two STNOs fully connected along their sides (Fig. 
H^aV). In contrast to circular vortex-based STNOs used, e.g., 
in J9|], we have chosen the square-shaped nanoelements, be- 
cause in this case the dipolar coupling caused by surface mag- 
netic charges on closely placed parallel edges of two such 
nanoelements should be much stronger than the correspond- 
ing coupling of two circular elements (we note in passing that 
the patterning of square- shaped nanoelements should be also 
more convenient from the technological point of view). For all 
three systems we have studied their dynamic behavior for cur- 
rents through the 1st nanocontact in the region 0.5 < I\ < 10 
mA, whereby the second current strength was kept constant 
(h = 5 mA). Simulations for each value of I\ were performed 
independently. 

The choice of the current and spin polarization directions 
for the two contacts system requires a special discussion. The 
current direction defines the rotation sense of its Oersted field 
and thus - the rotation sense of the equilibrium Landau do- 
main pattern in each nanosquare. The latter, in turn, defines 
the chirality of the corresponding vortex (polarities of all vor- 
tices are identical due to the presence of the out-of-plane ex- 
ternal field). If the current directions are the same for both 
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FIG. 5. (color online) Coupling energy vs current strength for sys- 
tems of two vortex STNOs: triangles - system with dipolar coupling 
only (Fig. 01)), squares - 'bridge' configuration (Fig. Ha)), circles 
- full exchange and dipolar coupling (Fig. |Ha)). 



contacts, then the sign of the out-of-plane spin polarization 
projections P pQrp is also the same (we remind, in a real sys- 
tem the sign of P perp is determined by the current direction 
and the sign of the out-of-plane magnetization component of 
the 'hard' magnetic underlay er). In this situation the gyration 
sense of both vortices coincide, providing favorable condi- 
tions for the STNOs synchronization. However, if the rotation 
senses of equilibrium domain patterns are the same in both 
nanosquares, an antivortex is formed between the vortices in 
exchange-coupled systems (cases (ii) and (iii) listed above). 
Such an antivortex is a strong undesired perturbing factor, 
because its dynamics can not be satisfactorily controlled by 
external parameters even in the regime where no additional 
V-AV pairs are generated: the antivortex annihilates unavoid- 
ably with one of the oscillating vortices (after the V-AV dis- 
tance reaches some minimal value), thus destroying one of 
STNOs. For these reasons we propose to use currents having 
opposite directions for adjacent point contacts, so that the ro- 
tation senses of the vortices created by corresponding Oersted 
fields are also opposite; for this configuration no antivortex 
is formed in exchange-coupled systems. In this case the spin 
polarization directions would also be opposite, so that the gy- 
ration senses of adjacent vortices remain the same. 

Results for the case (i) of dipolar-coupled nanoelements are 
summarized in Fig. [2] For the current directions and spin 
polarizations used here, in the steady precession regime both 
vortices rotate clockwise around the contacts. After the initial 
in-plane field pulse expels the vortices out of contact areas, the 
magnetization within these areas rotates 'in-phase' for several 
tens of nanoseconds, but afterwards a transition to an 'out-of- 
phase' rotation occurs, as shown in Fig. 12b). Details of this 
highly interesting transient dynamics governed by interaction 
between vortices will be discussed elsewhere. 

From the oscillation power spectra of the sum m* ot (t) = 
m*(t) + m*(0 of magnetization components under the 1st and 



2nd contacts shown in Fig. EJc), it is clear that the synchro- 
nization of two STNOs occurs nearly in the entire current in- 
terval of the regular vortex precession (i.e., where no V-AV 
pairs are generated) 1.5 < I\ < 7 mA. In this interval, the to- 
tal power at the fundamental frequency / = 0.50 GHz nearly 
disappears (for 3 < I\ < 6.5 mA it is below 10" 6 of the power 
for a single STNO), thus signaling the perfect out-of-phase 
synchronization of two oscillators. Indeed, the snap-shot in 
Fig. |2a) shows that the vortex cores precess in-phase, but 
due to their opposite chiralities, magnetic moments under the 
contacts rotate out-of-phase, so that their sum m* t is approxi- 
mately zero. 

In an experimental system, in order to obtain the maximal 
output power, one would either perform the phase shift of one 
of the currents by 7r, or use the underlayers with the opposite 
orientations of the in-plane magnetization projections; in both 
approaches the GMR signals from the two contacts would be- 
come in-phase. In this case one would obtain the dependence 
of the total oscillation power on the current strength P(h) 
shown in Fig. Etd), where we have plotted the total oscillation 
power of m x m (t) = m\(t) + m*(0 • exp(/7r) = m\(t)-m x 1 (t). As it 
should be, we observe a strong enhancement of P in the syn- 
chronized regime. This enhancement is even somewhat more 
than two times (the maximal increase expected for two ideal 
synchronized oscillators), because for low and high currents, 
where the oscillators are not synchronized, the magnetization 
under contacts is not entirely coherent. 

For the second system, where the nanosquares are con- 
nected via a bridge, vortex STNOs interact both via the dipolar 
and exchange interaction. We remind, that due to the chosen 
current directions we avoid the formation of the antivortex 
between the vortices, so that the we still can obtain a stable 
steady-state precession of vortices and their synchronization, 
as shown in Fig. [3a). Due to the presence of an additional 
exchange coupling, synchronization starts to establish itself 
already for the smallest current where the vortex precession 
could be induced (see Fig. [3l (b) and (c)). As in the previ- 
ous case, magnetization oscillations under the contacts in the 
synchronized regime are out of-phase. 

The third system - fully connected square-shaped nanoele- 
ments (so that they form a single rectangle) - demonstrates a 
qualitatively different type of synchronization (Fig. ffla): the 
strong exchange interaction between STNOs leads mostly to 
out-of-phase vortex core oscillations (and the vortex trajec- 
tories are strongly non-circular), so that magnetization oscil- 
lations under the contacts are in-phase. However, the syn- 
chronized state observed on previous systems - with the in- 
phase vortex oscillations and out-of-phase m x -oscillations un- 
der the contacts - is also possible, as it can be seen from the 
spectral power 'gaps' around I\ = 3.5 mA and 6.5 mA in 
Fig. Sfb). The existence of such multiply synchronized states 
has been reported for another system of synchronized STNOs 
W22I1 . We also point out that for large currents, where in previ- 
ous systems the V-AV pairs were generated, we still observe 
here synchronization, but with several frequency jumps (Fig. 
Ub)); corresponding processes will be discussed in details 
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FIG. 6. Snapshots of the in-plane magnetization distributions for 10 nanooscillators in the transient (t = 15 ns, upper panel) and synchronized 
(t = 1 16 ns, lower panel) regimes. 
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FIG. 7. Total oscillation power (normalized to the oscillation power 
of one STNO) of the alternating sum of the in-plane components 
inside the point contacts m x ^{t) - 2;(-l) z_1 m*(0, which takes into 
account the opposite oscillation phases of the magnetization in adja- 
cent contacts. 



elsewhere. 

One of the most important parameters of a synchronized 
system is the coupling energy between the oscillators E coup . 
This quantity defines not only the current range where syn- 
chronization occurs and the maximal allowed frequency mis- 
match between single oscillators |4|,|6|,|9|], but also the stability 
of the synchronization against thermal fluctuations. For some 
simple cases there exist analytical and semianalytical meth- 
ods for the estimation of E conp |4j,|9|], but for complex systems 
like dipolar- and exchange-coupled vortices reliable analytical 
methods to compute this quantity are lacking. 

For this reason we suggest here a purely numerical and gen- 
erally applicable method for calculation of the coupling en- 
ergy. First we note, that the total energy of the synchronized 
system even in a steady-state oscillation regime is not con- 
stant, but oscillates slightly in time: E tot = E tot (t). Averaging 
£tot(0 over at least one oscillation period (better over several 
periods), we obtain the averaged total system energy (E tot ). 
Subtracting from this quantity the equilibrium system energy 
E Qq at T = (where the Oersted field is already present, but 
the spin-transfer torque is still switched off), we obtain the 



dynamical contribution to the total energy E S y^ c , pumped into 

the system by the spin-transfer torque: E^ c = (^tot) _ E Qq . 
Then we repeat this procedure for each oscillator indepen- 
dently, obtaining the dynamical additives E^ yn and E^ yn for 
the oscillators 1 and 2 operating separately; all these ener- 
gies depend, of course, on the current strength. The difference 
between the dynamical additive for the system of vortices in 
the synchronized regime E^ c an d the sum of E dyn of two 
STNOs separately (operating at the same current strengths as 
the synchronized system) provides a dynamical coupling en- 
ergy: £"coup — E S y nc — {E y + E^ ). 

Results of this calculation for all three systems studied 
above are presented in Fig. [5] where the plots of the abso- 
lute values of £ C oup (which are by itself obviously negative) 
are shown; note the logarithmic scale of the energy axis. The 
numbers on the plot are the coupling energies for l\ = 5 mA 
(h - 5 mA, as always). First we note, that even for the system 
with only the dipolar interaction between STNOs the coupling 
energy is roughly twenty times larger than kT. For this reason 
we expect the effect of thermal fluctuations in these systems to 
be irrelevant. We have verified this expectation by micromag- 
netic modeling including thermal fluctuations at T = 300 K. 
The spectral amplitude map P(f, I\) for this temperature was 
very similar to that presented in Fig. |2c) for T = 0. Compar- 
ison of the total powers for T = and T = 300 K (Fig. [2d)) 
demonstrates that thermal fluctuations do not really affect the 
synchronization dynamics even for this system with the weak- 
est coupling between STNOs. We also note, that the coupling 
energy E coup = 8.4 x 10" 13 for this our system is roughly four 
times larger than ^coup = 2.3 x 10" 13 erg, calculated in 09D 
for two vortex-based STNOs with the circular shape; this dif- 
ference, however, should be attributed partially to the larger 
thickness of our nanoelements. 

For the fully exchange-coupled system the relation of typ- 
ical coupling energy to the thermal energy is E coup /kT « 
2.9 x 10- n /4.1 x 10" 14 « 700 (Fig. 0), so that in this sys- 
tem thermal fluctuations should play no role whatsoever. 

The last very important point is that our design allows 
the synchronization of an arbitrary number of vortex-based 
STNOs, arranged, e.g., as a regular chain. To demonstrate this 
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fact, we have performed simulations of a system consisting of 
10 dipolarly coupled oscillators placed as shown in Fig. the 
total system size is 5430 x 500 nm 2 . Here, total currents (not 
the current densities) flowing through all contacts have been 
taken to be equal; the value of this current was varied from 0.5 
to 7 mA. 

In these simulations, we have addressed another very im- 
portant problem: how stable is the synchronization with re- 
spect to random fluctuations of the point contact diameters, 
unavoidable in technological processes used for manufactur- 
ing of such contacts. Simulating a single STNO, we have 
found that the vortex gyration frequency is nearly indepen- 
dent on the point contact diameter, if the total current through 
the contact is kept constant. The reason is, that the vortex os- 
cillation orbit lies outside the point contact region, where the 
Oersted field (and thus - the vortex dynamics) depends only 
on total current flowing through the contact. Hence a robust 
synchronization of many STNOs based on point contacts with 
very different diameters can be expected, as long as the cur- 
rents strengths (which can be controlled experimentally very 
precisely) through all contacts are the same. In our simula- 
tions we have used a system of 10 STNOs, with contacts di- 
ameters having a normal distribution with the mean D = 90 
nm and dispersion cr = 30 nm; in our concrete system realiza- 
tion the contacts have the diameters 60 < D < 120 nm. Fig. 
[6] displays in-plane magnetization distributions for two time 
moments: a few nanoseconds after the external field pulse has 
been switched off it = 15 ns) and 100 ns later, when oscillators 
are almost perfectly synchronized. Plot of the total oscillation 
power (normalized to the oscillation power of one STNO) vs 
current (Fig. [7]) shows that a huge power gain of nearly 100 
times (close to a maximal theoretical value) can be obtained 
in this highly non-ideal system for currents Iq > 4 mA. 

In summary, we predict, that in various kinds of systems 
consisting of several vortex-based STNOs driven by a spin- 
polarized current which flows through point contacts attached 
to square- shaped magnetic nanoelements, a highly stable syn- 
chronization of these nanooscillators can be achieved. Using 
a general method to calculate the coupling energy E coup of 
synchronized STNOs suggested by us, we show that in such 
a system E coup can be up to three orders of magnitude higher 
than kT, so that the influence of thermal fluctuations is ex- 
pected to be fully negligible. Finally, we demonstrate, that a 
robust synchronization of an arbitrary number of spin-torque 
driven vortex nanooscillators can be easily obtained using our 
suggested system setup. This synchronization remains stable 
also in a system where a significant (experimentally unavoid- 
able) distribution of the point contact diameters is present. 
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